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Abstract 

A new molecular motor is conceptually constructed from a double-walled carbon nanotube 
(DWNT) consisting of a long inner single-walled carbon nanotube (SWNT) and a short outer 
SWNT with different chirality. The interaction between inner and outer tubes is the sum of the 
Lennard-Jones potentials between carbon atoms in inner tube and those in outer one. Within 
the framework of Smoluchowski-Feynman ratchet, it is theoretically shown that this system in an 
isothermal bath will exhibit a unidirectional rotation in the presence of a varying axial electrical 
voltage. Moreover, the possibility to manufacture this electrical motor from DWNT is discussed 
under the current conditions of experimental technique. 

PACS numbers: 61.46,+w, 85.35.Kt 
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In order to illustrate that the second law of thermodynamics cannot be violated, Feyn- 
man introduced an imaginary ratchet-and-paw system in his famous lectures.— This system 
is now called Smoluchowski-Feynman ratchet because it was first described by Smoluchowski 
in 1912. 2 One can observe the ratchet effect that the ratchet rotates unidirectionally if two 
conditions are satisfied 3 -: (i) breaking of spatial inversion symmetry and (ii) breaking of 
thermal equilibrium. Based on Smoluchowski-Feynman ratchet, many statistical models 
are put forward, such as on-off ratchets,^ fluctuating potential ratchets,— fluctuating force 
ratchets,— £ temperature ratchets 9 - ' 10 ^ 1 and so on. These models can explain the directional 
matter transport at the molecular levels in biological systems with the aid of the special pro- 
tein systems . 12 i 13 i 14 Now the small devices or systems at molecular levels are called molecular 
motors if they can transform chemical, electrical or other forms of energy into mechanical 
energy. 

In another classic talk given by Feynman at the annual meeting of the American Phys- 
ical Society in 1959, he predicted that Human would manufacture devices even at the 
molecular levels in the future.-^ 5 - His prediction are being realized with the development 
of nanotechnology.— Especially, after the discovery of carbon nanotubes,— researchers de- 
signed many nanodevices, such as nanotube gears , 1§l19 nanotube oscillators , 20 i 21 i 22 i 23 nan- 
otube drills 24,25 and so on. Tuzun predicted a doped nanotube motor driven by laser,— 
and one of the present authors and Ou-Yang proposed a molecular motor constructed from 
a double-walled carbon nanotube (DWNT) driven by temperature variation.™ Kang et al. 
confirmed these predictions by molecular dynamics simulations and then put forward nan- 
otube motors driven by fluid.— Viewed from statistical mechanics, the DWNT motor in 
Ref. is a kind of temperature ratchet o 9llQl11 in essence. But it is not easy to control the 
variation of temperature. Thus we may ask: can we design a more convenient way, for 
example using the electrical field, to drive this motor? 

In this report, we show it is possible to construct a DWNT motor driven by electrical 
field in physical point of view. The motor consists of a long inner tube and a short outer 
tube. For simplicity, the interaction between inner and outer tubes is taken as the sum 
of the Lennard- Jones potentials between carbon atoms in inner tube and those in outer 
one. Within the framework of Smoluchowski-Feynman ratchet, we theoretically reveal that 
this system in an isothermal bath will exhibit a unidirectional rotation in the presence of a 
varying axial voltage. Moreover, we also discuss the possibility to manufacture this electrical 
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motor under the current conditions of experimental technique. 

As illustrated in Fig. Q the DWNT consists of a long inner tube and a short outer tube 
with different chirality. Without losing generality, we take (8, 4) and (14, 8) SWNTs, and 
put two-unit cells for outer tube. The SWNTs have the same axis. One end of inner tube is 
fixed while another is simply sustained. We put the DWNT in between, without contacting 
to, a pair of parallel electrodes on which we apply a voltage varying with time. The axis of 
DWNT is perpendicular to the electrodes. The whole system is put in a thermal bath, e.g., 
helium gas. Because of the extremely high Young's modulus 29 and stereo effect of carbon 
nanotubes, the collisions of helium atoms at low temperature, e.g., T = 50 K almost have 
not effect on the shape of carbon nanotubes as well as the coaxiality of two nanotubes. Two 
degrees of freedom will be excited by the collisions: the rotation of outer tube around the 
inner one and the slide of outer tube along the axis. In following contents, we assume that 
some device can fix the sliding degree of freedom of outer tube. The interactions between 
atoms in outer tube and those in inner tube provide potentials asymmetric with respect 
to the relative rotation of the two nanotubes because of the difference in chirality. This 
property breaks the spatial inversion symmetry and naturally makes our device satisfy the 
first condition to exhibit the ratchet effect. 

The applied varying electrical voltage will achieve the second condition of ratchet effect 
as discussed below. Recently Guo et al. found a large axial electrostrictive deformation 
of SWNTs in the presence of axially electrical field by using numerical Hartree-Fock and 
density functional calculations^ although SWNTs are usually regarded as non-piezoelectric 
materials. According to these authors, the electrical field changes the electronic structures of 
SWNTs and induces elongations of carbon-carbon bonds without changing bond angles. 39 
Although their calculations are concerned to the electrostatic field, it is easy to see that 
their results are also available to the varying electrical field with period in the order of 
nanoseconds because the response time of electron to electrical field is usually about several 
femtoseconds. Thus, bond elongations are synchronous with the variation of electrical field, 
which induces the shaking interaction between outer and inner tubes. Consequently, we 
obtain a fluctuating potential which, according to Ref. is sufficient to clear the second 
condition for ratchet effect in the statistical point of view. 

To specify the above idea, we select an orthogonal coordinate system Oxyz by adopting 



the convention in Ref. 



31 



whose 2-axis is the tube axis parallel to the translation vector^ of 
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SWNT, and x-axis passes through some carbon atom in inner tube. The angle rotated by 
outer tube around inner tube is denoted by 9. All coordinates of atoms are also calculated by 



the method in Ref. 
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For theoretical simplicity, the interaction between any atom % in outer 
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tube and atom j in inner tube is taken as Lennard- Jones potential^ u(rij) = 4e[(a/r, 
(a/rij) 6 ], where is the distance between atoms % and j, e = 28 K, and a = 3.4A.— 
The inner tube is long enough to be thought of as an infinite one. The total interaction, 
V(9, a cc ) = u ( r ij)i between outer tube and inner tube can be calculated by the method in 
Ref. l3lL where a cc is the length of carbon-carbon bond whose value is 1.42 A for zero voltage. 
Obviously, the potential can be expanded by Taylor series V(9,a cc ) ~ Vo(0)[l + a(6)e] for 
small elongations of bond-length, where V$(0) = V(9, 1.42A), a (9) is a function of 9 and e 
the bond elongation ratio. In Fig. |2l we show numerical results of the interaction between 
outer tube and inner tube for different bond-lengths. Here the potentials V(0, a cc ) are put 
to zero by subtracting some constants which have no effect to our final results. We find that 
the potentials are functions with period 7r/2 and that oi{9) is so weakly dependent on 9 that 
we can fit it by a(9) = —14.2 from Fig.|21 On the other hand, from Fig. 2 in the work of Guo 
et al., we know the bond elongation ratio shows approximately a linear dependence on the 
magnitude of electrical field with the slope (3 = 0.25 A/V provided that the applied voltage 
field is not too high. Moreover, the field dependence of bond elongation ratio is insensitive 
to the diameters and chiral angles of SWNTs. Therefore, we obtain a shaking interaction, 
V(9, t) = Vq{9){1 + /i[l + cos(27rt/T)]}, between outer tube and inner tube under a varying 
voltage U(t) = U [l + cos(2irb/T)] between two electrodes, where U , T and \i = aj3U /l are 
constant quantities. Here / is the distance between two electrodes and T is assumed greater 
than 0.1 ns. 

Next we will discuss the motion of outer tube in helium gas. Because the mass of outer 
tube is much larger than that of helium atom, its motion can be described by the Langevin 
equation mR 2 9 = —V'(9,t) — r]9 + where m ~ 10 -23 Kg, the mass of outer tube 

containing 496 carbon atoms, and R = 7.75A is the radius of outer tube. r\ is the rotating 
viscosity coefficient, and the dot and the prime indicate, respectively, differentiations with 
respect to time t and angle 9. £(£) is thermal noise which satisfies (£(£)) = and the 
fluctuation-dissipation relation (£(£)£(0)) = 2i]T5(t), u where T = 50 K and the Boltzmann 
factor is set to 1. We estimate mR 2 9/(i]9) < 10 -3 by taking 9 ~ 9/T and T > 0.1 ns. Thus 
it is reasonable to neglect the inertial term mR 2 9. In this case, the Fokker-Planck equation 
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corresponding to the Langevin equation is expressed as^ 
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where P(9, t) represents the probability of finding the outer tube at angle 9 and time t which 
satisfies P(9 + n/2,t) = P{6,t) and ft /2 P(6,t)d8 = I. It follows that the average angular 
velocity of outer tube in the long-time limits 
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Let D = — , t = D 1 t, T = D 1 J', Vq{6) = u(9)T, we obtain the dimensionless forms of 
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Because V (9) is a periodic function, we can expand it by Fourier series and find that it is 
well fit by V (9) = 29.12-0. 18cos40-16.14cos 80-12.53 cos 120-4.48 sin 45-19. 52sin 89 + 
11.52 sin 120 (K), where high-order terms are neglected because their coefficients are very 
small. Thus u\9) = V£(0)/T = T l l=i( v ke 4ike + vle~ 4ike ) with Vl = -0.179 - 0.007*, v 2 = 
— 1.562 — 1.291? and V3 = 1.382 — 1.504z. In the the long-time limit, P(9,t) can be expanded 
by Fourier series P(9,r) = X^ m =-ooP fimex P^[(2 n7rT 7 + 4m0]. Substituting them into 
(J3J), we arrive at a recursion equation 
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Considering the constraint p nQ = (2/n)5 n o coming from J^ 2 P(9,t)d9 = 1, we solve 
the recursion equation (JSJ) with [i = —0.01 for different dimensionless period J7, and then 
calculate the dimensionless average angular velocity by ©. The dimensionless period 
dependence of average angular velocity of outer tube rotating around inner tube is shown 



in Fig. 01 from which we find that: (i) — > if J — > oo which is equivalent to the 
case of an electrostatic voltage; (ii) (^) = 0if l 7 = which corresponds to the case that 
the voltage varies so quickly that outer tube can not respond to it in time; (iii) Outer tube 
rotates counterclockwise around the tube axis for some J , which corresponds to the positive 
value of (^), and clockwise for other J which corresponds to the minus value of (fj:); (iv) 
{^) has a significant value —347 mad when J c = 0.15. 

For helium at T = 50 K, we can calculate rj = 1722 K-ns from its value at 273 K.— It is 
necessary to emphasize that 77 is the rotating viscosity coefficient which is different from the 
viscosity coefficient of gas r] c in the common sense. The simple relationship between them 
is rj = 2irR 2 Lri c for cylinders, where L = 27.39 A is the length of outer tube. Therefore we 
obtain (0) = —10.08 nrad/ns (i.e., about one and a half rounds per second) when the period 
T c = 5.17 ns. Here we obtain these values in disregard of the fiction between nanotubes 
because it is very small.— 

The above discussions demonstrate that we can construct an electrical motor from DWNT 
in principle. Two key points are: (i) DWNT with different outer and inner chirality induces 
a potential that breaks spatial inversion symmetry; (ii) Some mechanism (varying voltage in 
this paper) induces the shaking of the potential that breaks thermal equilibrium. Therefore, 
similar results are available in any DWNT with different outer and inner chirality although 
our discussions are focused on (8,4)@Q4,8) tube. Additionally, we use the results obtained 
from density functional theory in Ref. ; 30 that shows large axial electrostrictive deformation 
for armchair and zigzag tubes with different diamters. These investigated cases reveal that 
the electrostrictive deformation is insensitive to the diameters and chiral angles of SWNTs. 
Generally speaking, different voltage dependence of bond-elongation might exist for different 
chirality. Even if were the case, it does not hinder from inducing the shaking of the potential. 
In short, different DWNT and different voltage dependence of bond-elongation do not change 
the qualitative conclusion of this paper but merely modify the quantitive value of the average 
angular velocity. 

Can this device be made in reality? There are several technological difficulties. The first 
one is to make the DWNT with different chirality of outer and inner tubes. Experimentalists 
can overcome it by producing DWNT based on the synthesis technique in Ref. |3jJ The 
second difficulty is the varying voltage. Our device requires voltage (Uo ~ 3 V for / ~ 0.1 
/iin) with high frequency of gigahertz. The third one is the method to restrict the sliding 
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degree of freedom of outer tube. It is possible to manufacture the DWNT motor driven by 
varying electrical voltage provided that researchers overcome the above difficulties. 

In summary, we conceptually design a new molecular motor from DWNT driven by 
varying electrical voltage. In this sense, it can be called "electrical motor" in nanometer 
scale. Similarly, we may design another molecular motor if we replace carbon nanotubes 
with boron nitride nanotubes because the latter is a piezoelectric material.— People have a 
long dream to make machines in nanoscales. The key component to these machines is the 
power device. Our motor is expected to be the power device if it is manufactured in the 
future with the development of nanotechnology. 

We are grateful to Dr. Q. X. Li and B. Y. Zhu for their technical helps. 
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FIG. 1: Schematic diagram of the DWNT motor driven by varying electrical field. The DWNT 
consists of a long inner tube and a short outer tube. It is put between a pair of electrodes with a 
varying voltage U (t) applied on. 
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FIG. 2: Interactions between outer tube and inner tube with different lengths of carbon-carbon 
bonds. Here the potentials V(0, a cc ) are put to zero by subtracting some constants. 
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FIG. 3: Average dimensionless angular velocity (dO/dr) of outer tube rotating around inner one in 
isothermal bath when varying axial voltage with dimensionless period J is applied. The minus sign 
means the clockwise rotation around z-axis while the positive represents counterclockwise rotation. 
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